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Abstract 

We present a versatile density functional approach (DFT) for calculating the 
depletion potential in general fluid mixtures. In contrast to brute force DFT, 
our approach requires only the equilibrium density profile of the small parti- 
cles before the big (test) particle is inserted. For a big particle near a planar 
wall or a cylinder or another fixed big particle the relevant density profiles are 
functions of a single variable, which avoids the numerical complications inher- 
ent in brute force DFT. We implement our approach for additive hard-sphere 
mixtures, comparing our results with computer simulations for the depletion 
potential of a big sphere of radius Rb in a sea of small spheres of radius Rs 
near i) a planar hard wall and ii) another big sphere. In both cases our results 
are accurate for size ratios s = Rs/Rt as small as 0.1 and for packing fractions 
of the small spheres rjs as large as 0.3; these are the most extreme situations 
for which reliable simulation data are currently available. Our approach sat- 
isfies several consistency requirements and the resulting depletion potentials 
incorporate the correct damped oscillatory decay at large separations of the 
big particles or of the big particle and the wall. By investigating the deple- 
tion potential for high size asymmetries we assess the regime of validity of 



the well-known Derjaguin approximation for hard-sphere mixtures and argue 
that this fails, even for very small size ratios s, for all but the smallest values 
of r]s where it reduces to the Asakura-Oosawa potential. We provide an ac- 
curate parametrization of the depletion potential in hard-sphere fluids which 
should be useful for effective Hamiltonian studies of phase behavior and col- 
loid structure. Our results for the depletion potential in a binary hard-sphere 
mixture, with size ratio s = 0.0755 chosen to mimic a recent experiment on 
a colloid-colloid mixture, are compared with the experimental data. There is 
good overall agreement, in particular for the form of the oscillations, except 

at = 0.42, the highest value of packing fraction considered. 
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I. INTRODUCTION 



Two big colloidal particles immersed in a fluid of smaller colloidal particles or non- 
adsorbing polymers or micelles experience an attractive depletion force when the separation 
h of the surfaces of the big particles is less than the diameter of the small ones. The expulsion 
or depletion of the small particles gives rise to anisotropy of the local pressure which results in 
the effective attractive force between the big ones. Asakura and Oosawa and, independently, 
Vrij used excluded volume arguments to determine the effective potential between two big 
hard spheres (modeling the colloids) assuming that the small particles or polymers form a 
mutually non-interacting fluid whose centers are excluded from the surfaces of the colloids 
by a distance Rs The resulting depletion potential is attractive for h < 2Rs and is 
zero for h > 2Rs\ it increases monotonically with h from its value at contact, /i = 0, 
and is proportional to r/^, the packing fraction of the small particles [see, c.f., Eq. (0)]. 
Much attention has been paid to depletion induced attraction within colloid science since 
it provides an important driving force for phase separation and fiocculation phenomena in 
mixtures of colloids and in colloid-polymer mixtures. From a statistical mechanics viewpoint 
depletion forces are of considerable interest since they arise from purely entropic effects 
because the bare interactions between the particles are hard-sphere-like. Formally, it is 
the integrating out of the degrees of freedom of the small particles which gives rise to the 
effective interaction between two big ones. 

Although colloid-polymer mixtures can, under favorable circumstances, be modeled by 
a binary mixture of hard spheres and ideal, non-interacting polymers (we term this the 
Asakura-Oosawa model), for mixtures of colloids or colloids and micelles a more appropriate 
zeroth-order model is a binary hard-sphere mixture, i.e., the small particles are not inter- 
penetrating but experience mutual hard-sphere repulsion. In this case it becomes a key 
question as to how the depletion potential between two big hard spheres is influenced by 
interactions between the small spheres. For high packing fractions r^^, one might suppose 
that the small spheres exhibit pronounced short-ranged correlations (layering) leading to 
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significant changes in the depletion potential. This would, in turn, have repercussions for 
the phase behavior of the bulk mixture, making this significantly different from that of 
the Asakura-Oosawa model. Such considerations have prompted several recent theoretical 
investigations of phase behavior based on an effective one-component depletion potential 
description of model colloidal mixtures [0-0]. The crucial ingredient in such investigations 
is an accurate depletion potential. 

Having a proper understanding of depletion potentials is not only relevant for bulk phase 
behavior; it is of intrinsic interest. In recent years a variety of experimental techniques have 
been developed which measure, directly or indirectly, the depletion potential between a 
colloidal particle, immersed in a sea of small colloids or polymers, and a fixed object such 
as a planar wall ||^. Video microscopy has also been used to determine depletion forces for 
a single big colloid in a solution of small colloids inside a vesicle - a system which resembles 
hard spheres inside a hard cavity 0. Very recently Crocker et al. measured the depletion 
potential between two big PMMA spheres immersed in a sea of small polystyrene spheres 
for a range of packing fractions of the latter (see, c.f. Sec |V B|) . At low values of 1]^ the 



measured depletion potential is well-described by the Asakura-Oosawa result but at higher 
packing fractions the potential exhibits a repulsive barrier and for rjs > 0.26 the depletion 
potential is damped oscillatory with a wavelength that is of the order of the small particle 
diameter. As experiments grow in sophistication and in resolution it is likely that further 
details of depletion potentials will be revealed whose interpretation will require a reliable 
and versatile theoretical approach. Such an approach should be able to tackle experimental 
situations where t]s is rather high and to treat general 'confining' geometries. The latter 
include a big particle near a planar wall or in a wedge or cavity, as well as the case of a 
big particle near another, fixed big particle. In this paper we describe such a theory for the 
depletion potential based on a density functional treatment (DFT) of a fluid mixture. Our 
treatment avoids the limitations of the virial expansion (in powers of rjs) and the uncontrolled 
nature of the Derjaguin approximation which are inherent in recent approaches to 
depletion forces in hard-sphere mixtures. It is less cumbersome than the alternative integral 



equation treatments |TO|,|TT| and more easily adopted to different geometries. A l<;ey feature 
of our treatment is that it does not require tlie calculation of the total free energy of the 
inhomogeneous fluid or of the local density of the small particles in contact with the big 



particle [0,|T^]. The method is much easier to implement than a direct minimization of 
the free energy functional which is numerically very demanding when any symmetry of the 
density profile of the small spheres is broken by the presence of the big particle. 



The paper is arranged as follows: Subsection |Il7\| defines the depletion potential in 
an arbitrary mixture, showing how this is related to the one-body inhomogeneous direct 
correlation function of the big particles. In Subsec. |IIB| we use this result to derive an 



explicit formula for the the depletion potential in the low-density limit, where the densities 
of all species approach zero. For the particular case of a binary hard-sphere mixture in 
this limit we recover the Asakura-Oosawa result. Subsection |II C| describes the general 



asymptotic behavior of the depletion potential for h oo while Subsec. [II D| describes the 
implementation of the theory for a binary hard-sphere mixture using the DFT of Rosenfeld 
| |13| . In Section |T| we present several comparisons of our hard-sphere DFT results, for 
both sphere-sphere and (planar) wall-sphere depletion potentials, with those of computer 
simulation [|I^,|r^. Our theory performs well for all size ratios s = Rg/Rb and packing 
fractions t]s for which simulation results are available. We show that the leading-order 
asymptotic result for the depletion potential provides an excellent account of the oscillations 
in the calculated potential not only at longest range but also at intermediate separations 
of the big spheres. Section |^ is concerned with assessing the regime of validity of the 
well-known Derjaguin approximation which relates the force between two big objects to the 
integral of the solvation force, or excess pressure, of the small particles confined between two 
planar walls [see, c.f., Eq. (|35|)]. We argue that this approximation is not reliable for the 
hard-sphere mixture even when the size ratio s is very small. In Subsec. |V A| we describe 



a simple but accurate parametrization of the depletion potential suitable for a big hard 
sphere near a planar hard wall and for the potential between two big hard spheres. Such a 
parametrized form should prove useful for effective Hamiltonian studies of phase behavior 



and colloid structure [g]-^ . Subsection [VB| presents results for the depletion potential in a 
binary hard-sphere mixture where the size ratio is chosen to mimic the system considered 
in the experiments of Ref. |^. We conclude in Sec.|VI| with a discussion and summary of our 
results. 



II. THE DEPLETION POTENTIAL 



A. General theory 



We consider a general mixture of v components in which each species i ii = 1, . . . , z/), 
characterized by its radius Ri, is coupled to a reservoir with chemical potential /ij and is 
subject to an external potential Vi{r). The mixture at thermodynamic equilibrium can be 
described by the set of number density profiles {pj(r)}. For such a mixture we wish to 
calculate the depletion potential, or the depletion force, between an object fixed at position 
ri and a second one fixed at r2. Without loss of generahty the position ri of the first object 
is chosen as the origin of the coordinate system. This fixed object then exerts an external 
potential on the particles constituting the mixture. The external potential can represent a 



planar hard wall [T^l or a fixed particle of the mixture, or more generally, a curved surface 
ir?! or soft planar walls [|TS[. If the depletion potential between two particles of the mixture 
is to be calculated either particle can be chosen to act as the external potential and this 
point will be addressed in more detail in a later section. 

In the following the second object is a test particle of a species denoted as b. The grand 
potential of the mixture when the test particle is fixed at the position r;, in the presence of 
the fixed object exerting the external potential ^(r) is denoted by fltbi't^b] {/^i};^)- Wt{rb), 
the quantity of interest here, is defined as the difference of grand potential between a con- 
figuration in which the test particle is in the vicinity of the fixed object and one in which 
the test particle is deep in the bulk, i.e., oo: 

Wtivb) = ^tb{TCb] {yUi}; T) - Vltbiri, oo; {/i^}; T). (1) 
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In order to calculate this difference the test particle can be moved along any path from 
one configuration to the other. A particular path which simplifies the calculation is via the 
reservoir. This path can be divided into two steps. In the first step the test particle is 
removed from the bulk at ^ oo and put into the reservoir. In the second step the test 
particle is taken from the reservoir and is inserted back into the mixture but now at r^. The 
formal means to describe particle insertion in a general mixture is the potential distribution 



theorem and we employ this in the grand ensemble [|1^ . 

The potential distribution theorem provides an expression for the partition function 
Stbirb, {/ii};T) of the mixture after a test particle of species b is inserted at position r;, in 
terms of the partition function of the mixture S({/ij};T) and the number density profile 
Pb(r) of species b before the particle insertion: 

Etbin; T) = exp {(SiVM - /x,)) ptin) T), (2) 

with f]~^ = ksT and A;, the thermal wavelength of species b. Together with a well-known 
result from density functional theory (DFT) 



Al p,(r) = exp - H(r)) + 4^^; M)) , (3) 

it follows that the one-body direct correlation function c[^'' of species b can be written as 

c'b\n] {pi}) = In (^tbin; {pi}; T)/H({/ii}; T) j 

= pn{{fi,}-T) ~ /3ntb{n;{fii};T), (4) 

i.e., —Pci^\rb] {pi}) describes the change in the grand potential of the whole system due to 
insertion of a test particle. The grand potential difference defined by Eq. (|I|) can now be 
expressed in terms of the difference of one-body direct correlation functions: 

PWM = c^^\n ^ oo; M) - ci;\n; M). (5) 

As the potential distribution theorem, Eq. (j^), is a general result, valid for any number 
of components, for arbitrary densities of all components and, in fact, for any inter-particle 



potential function the same generality holds for Eq. (|^). No approximations have been made 
so far. However, in order to use Eq. to calculate l3Wt{r) an explicit procedure that can 
treat a mixture must be applied. Simulations provide such a procedure as does density 
functional theory. We shall consider both here. 

We emphasize that the direct correlation function entering Eq. depends on the equi- 
librium density profiles before the test particle of species b is inserted at position r?,. This 
observation simplifies the calculation of l3Wt{r) dramatically because the symmetry of the 
relevant density profiles {pj(r)} is determined solely by the symmetry of the external poten- 
tials and therefore depends only on the nature of the object that is fixed at the origin. If this 
object is a structureless planar wall and in the absence of spontaneous symmetry breaking 
such as prefreezing or crystalline layer formation the density profiles of all species reduce to 
one-dimensional profiles {pi{z)} with z the distance perpendicular to the wall. For a fixed 
spherical or cylindrical wall or particle the density profiles {pi{r)} depend only on the radial 
distance. Even if the fixed object is a wall of more general shape, so that there is no simple 
symmetry involved in the problem, calculating the density profiles before particle insertion 
is much easier than after insertion, when the broken symmetry due to the presence of the 
test particle leads to a more complex dependence of the profiles on the coordinates. 

While Eq. can be evaluated for arbitrary densities of species b within the present 
DFT approach, a particular limit in which the density of species b goes to zero is considered 
now. This dilute limit is especially important since it arises in the context of measuring 
depletion forces and in formal procedures for deriving effective Hamiltonians for big particles 
by integrating out the degrees of freedom of the small particles. For example, if in a binary 
mixture the degrees of freedom of the small particles are integrated out the resulting effective 
one-component fluid can be described by an effective Hamiltonian containing a volume term, 
to which only the small particles contribute, a one-body term, in which a single big particle in 
a 'sea' of small particles contributes, a two-body term, a three-body term and so on . For 
highly asymmetric mixtures the most important contributions come from the volume and the 
one- and two-body terms. This assumption is substantiated by the results of calculations 
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of three-body contributions reported in Ref. [|T4| for a size ratio s = 0.1. Three-body 



contributions also seem to be small for s = 0.2 ||21|. Note that the two-body term describes 
an effective pairwise interaction potential between two big particles which turns out to be 
precisely the depletion potential, i.e., pWt{r) evaluated in the dilute limit 0. 

In the grand ensemble the dilute limit can be obtained by taking the limit in which the 
chemical potential of species b, fib —oo, with the chemical potentials of all other species 
{fiijLb} kept fixed. The depletion potential is then given by 

f3W{r) = lim (3Wt{r) 

which contains no explicit dependence on the external potentials that are present, i.e., the 
depletion potential depends only on the intrinsic change of the grand potential. 

Although in the dilute limit both the density profile Pb(r) and the bulk density pjj"''^ = 
Pfe(oo) of species b vanish, the ratio stays finite and the depletion potential can also be 
obtained from the result 

l3W{r) = - lim In ( -^^^ ) - H(r) + H(oo). (7) 



which takes a more familiar form if we re-write Eq. (|^ as p{r)/p{oo) = exp[—f3{W{r) + 
Vft(r))], where p(r) is the probability density of finding the particle of species 6 at a position 
r and we assume Vb{c)o) = 0. This route to the depletion potential was employed successfully 
in a grand canonical Monte Carlo simulation of a big sphere in a sea of small hard spheres 



near a hard wall [|T6|. It is also the route used to obtain W{r) from experiment |jT8|,| 
Note that in the same limit fib — ^ — oo the density profiles of all other species {pj^b(r)} 
reduce to those of a z/ — 1 component mixture. 

B. The low density limit 

In order to implement the formal result in Eq. a way of determining the direct 
correlation function cj^^ is required. It is convenient to adopt a DFT perspective. In density 



functional theory the intrinsic Helmholtz free energy functional can be divided into an ideal 
gas contribution plus an excess over the ideal gas contribution. While the former is known 



exactly, in general, only approximations are available for the excess part [|20[. One important 
exception is the excess free energy functional for a general mixture in the low density limit, 
i.e., in the limit of all densities going to zero. By means of a diagrammatic expansion it can 
be shown that the exact excess free energy functional in this limit is given by 

hm /?J-,,[{p,}] = -i^ / dV / dVp,,(r)p,(r')^(r-r'), (8) 

where fij is the Mayer bond between a particle of species i and one of species j. 

For a binary mixture in the low density limit the depletion potential acting on a big par- 
ticle b can be calculated from Eq. (|]) using the definition of the one-body direct correlation 
function given within density functional theory, 

and we obtain 

l3W{r) = - j dV (p,(r') - p.(oo))/,,(r - r'), (10) 

where s refers to the small particles. In the same limit the density profile of the small 
particles reduces to the density profile of an ideal gas in the external potential K(r), i.e., 
Ps(j) = Ps(oo) exp{—/3Vs{r)) and the depletion potential can be written as 

(3W{r) = -p,(oo) j dV (exp[-/3K(r')] - l)/fe.(r - r'). (11) 

This result is more familiar for the case of a binary hard-sphere mixture with sphere radii Rf, 
and Rs where /^^(r — r') = ~Q{{Rh + Rg) — |r — r'|), where G is the Heaviside function. Then 
Eq. (|l^) reduces to the well-known Asakura-Oosawa depletion potential As an example 
we consider the depletion potential between two big spheres; in this case exp[— /9V^(r)] — 1 = 
fbs{r) and the sphere-sphere depletion potential can be expressed as 

pwi'^ih) = -p,(oo)7r(2i?, - h) !^R,[R, + - ^[R, + ^] - for h< 2R, 

= for h> 2Rs, (12) 
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where h is the separation between the surfaces of the big hard spheres. As a second example 
we consider a big sphere near a planar, structureless hard wall. The depletion potential is 
then 

PW^Ii'W = -2ps{oo)n{2R, - h) ^^]- ^[R, - :|] - ^| for h < 2R, 

= for h> 2Rs, (13) 

where h is the separation between the surface of the big hard sphere and the hard wall. 

It is important to recognize that Eq. ([TT|) provides the exact low density expression for 
the depletion potential even if the interactions between the species or between the wall and 
species s are soft and possibly contain an attractive part so that the Mayer / functions 
cannot be expressed in terms of the Heaviside function G. In general there is also a direct 
interaction potential between two big particles, or between a single big particle and a wall, 
so that the total effective potential, after integrating out the degrees of freedom of the small 
particles, is the sum of the intrinsic contribution - the depletion potential - and the direct 
interaction potential Vb{r), i.e., 

<t>tot{r) = W{r) + V,{r). (14) 

For example the total effective potential between two big hard spheres in the sea of small 
hard spheres is ^tot{r) = W{r) + Vb{r), with Vb{r) the hard-sphere potential between the 
two big ones, and it is ^tot{f) which constitutes the effective pair potential in the effective 
one-component Hamiltonian for the big spheres 0. 

It is instructive to note that the functional given by the r.h.s. of Eq. (§) generates the 
appropriate depletion potential for the original Asakura-Oosawa model |]l| of a mixture of 
colloids and ideal, non-interacting polymers. This model binary mixture is specified by fee, 
fcp, and fpp, the Mayer / functions describing the pairwise interactions between two colloids, 
between a colloid and a polymer, and between two polymers, respectively, fpp is set to zero in 
order to describe the ideal, non-interacting polymer coils. The resulting depletion potential 
is still given by Eq. ([TT| ) but this result now holds for all polymer densities Ps(oo) not just 
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in the dilute limit Ps(c>o) — > 0, because the polymer is taken to be ideal. The total effective 
potential between two colloids is then given by Eq. (p!^), with the Asakura-Oosawa result 
(p!^) for W{r), which may be employed in an effective Hamiltonian for the colloids 0]. 

C. Asymptotic behavior 

In the previous subsection we showed that for the case of hard spheres the depletion 
potential reduces in the low density limit to the Asakura-Oosawa result. Examination of 



Eqs. (|T2D and ([131) shows that this potential is identically zero for separations h between the 
spheres or between the sphere and the wall that are greater than 2Rs = as- Outside the low 
density limit of the small particles this is no longer valid. From the general theory of the 
asymptotic decay of correlations p2[ it is known that for systems in which the interatomic 
forces are short-ranged, i.e., excluding power-law decay, the density profiles of both com- 
ponents of a binary mixture exhibit a common damped oscillatory form in the asymptotic 
regime, far from the wall or fixed particle, which is determined fully by the pole structure 
of the total pair correlation functions hij{r) of the bulk mixture. The depletion potential 
is related to the density profile of species b via Eq. and therefore its asymptotic behav- 
ior should be related directly to that of this density profile. In order to understand this 
connection in more detail we first recall some arguments from Ref. ||22|| . 



A bulk binary mixture consisting of small particles of density p^"''^ and big particles of 
density p^"''^ is considered. The total correlation functions in the bulk, hij{r), with i,j = s,b 
are related to the radial distribution functions gij{r) via hij{r) = gij{r) — 1 and to the two- 
body direct correlation functions c[^\r) via the Ornstein-Zernike relation for mixtures. In 
Fourier space the latter can be expressed as 

Uq) = (15) 
where hij{q) is the 3 dimensional Fourier transform of hij{r), the numerator is given by 
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NUq) = £i?(?)+pr(ci?(?)'-aiaHg)ai?(g)), 

NM = ci?(g) + p'f'ic^^Jiqr - cEHg)ai?(g)), (16) 



and 



D{q) = (1 - pT^'cfJiqm - Pr'4?(?)) - PrVr'ci?(?)^ (17) 

is a common denominator. The total correlation function in real space can be obtained by 
taking the inverse Fourier transform: 



oo 



rhiji^) ^ 2^ y '^^ gsin(gr) hij{q) (18) 







which can then be evaluated by means of the residue theorem. If g„ denotes the n-th pole 
in the upper complex q half-plane and i?„ the corresponding residue of qhij{q), the total 



correlation function can be written as 22 



rh.Ar) = ^Y.^''"'^^n. (19) 

n 

From this equation it becomes clear that the asymptotic behavior of hij (r) is dominated by 
the pole or poles qn with the smallest imaginary part, since this gives rise to the slowest 
exponential decay. 

For all pairs i,j = b,s the poles are determined by the condition D{q) = p2|. For a 



binary mixture in the dilute limit of the big particles, i.e., p^"''^ 0, the general theory of the 
asymptotic decay simplifies considerably and from Eq. ([T7| ) we see that the pole structure 
of all three total correlation functions can be obtained from the solutions of the equation 

l-pf'4?(9)=0, (20) 

with cis^(g) referring to the fluid of pure s at density p^"''^. In general there will be an infinite 
number of solutions of Eq. (^), but only the solution g„ = g = ai + iao with the smallest 
imaginary part is important for the following. The asymptotic behavior of the radial 
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distribution functions Pi{r)/pY^^ of species i around a fixed particle, which can be a small 
or a big one, can be ascertained and it follows that the density profiles exhibit asymptotic 
decay of the form 

Piir)-p^^^ ^exp(-aor)cos(air - 0pi), r ^ oo, (21) 

r 

with a common characteristic inverse decay length and wavelength of oscillations 27r/ai 
for both species i = s,b. Remarkably, exactly the same inverse decay length and wavelength 
also characterize the asymptotic decay of the density profiles close to a planar wall. This is 
given by ||,|3 



Pi(z) - p-"'^ ~ A^iexp(-ao2;) cos(ai2; - 6^i), ^ ^ oo, (22) 

for i = s,b. The amplitudes Api and A^i and the phases Bp, and Q^i do depend on species 
i and whether a particle or wall is the source of the external potential. Note that from 
Eq. (^) it follows that in the dilute limit for the big particles Oq and Oi are functions of the 
packing fraction of the small particles only; thus they do not depend on the size ratio. 

The asymptotic behavior of the depletion potential can now be obtained from Eq. (^. 
Assuming that the external potential acting on the big spheres is of finite range the depletion 
potential between two big spheres has an asymptotic behavior of the form 

pWir) ^ - In (^1 + ^ exp(-aor) cos(air - ep,)/p'^'^^ 

^_il!^exp(-aor)cos(alr-epfe)M"'^ r ^ oo, (23) 
r 

and that between a single big sphere and a planar wall takes the form 

pW{z) A^b exp{-aoz) cos{aiz - Q^b)/pT^, z ^ oo. (24) 

We shall see later that our DFT results for the density profiles and the depletion potential 
conform with these asymptotic results at very large separations and, strikingly, at interme- 
diate separations. 
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D. A density functional approach for hard spheres 



For the system of primary interest, namely the mixture of hard spheres, a very rehable 
DFT exists, namely the Rosenfeld fundamental measures functional [1^. While, in principle, 
this functional also can treat generally shaped convex hard particles [^], its application has 
been restricted to the particular cases of hard spheres and parallel hard cubes 

In the low density limit the Rosenfeld functional reduces to the exact excess free energy 
functional of Eq. (^). For arbitrary densities it has the following structure: 

= J dV $ (K(r)}), (25) 

where $ is a function of a set of weighted densities {ua} which are defined by 



i=l 



j d'r'p,{v')u;t{v-v'). (26) 



The weight functions cjf in Eq. (|26|) depend only on the geometrical features, the so-called 
fundamental measures, of species i. Explicit expressions for the weight functions of hard- 



sphere mixtures and for $ can be found in Ref. |13] and in Ref. [26|. The Rosenfeld functional 



has the following properties: i) the free-energy of the homogeneous mixture is identical to 
that from Percus-Yevick or scaled-particle theory and ii) the pair direct correlation functions 
of the homogeneous hard-sphere mixture, generated by functional differentiation of J^ex-, are 
identical to those of Percus-Yevick theory. The index a labels 4 scalar plus 2 vector weights 
r3| . While the on^fma/ functional given in Ref. |T3] did not account for the freezing transition 



of pure hard spheres, more sophisticated extensions do account for freezing; the weight 
functions remain the same but $ is changed slightly. For the depletion potential problems 
under consideration the different versions give almost identical results for bulk packing 
fractions rjs < 0.3 |2^. At higher packing fractions the density profiles of the small spheres 
Ps{z) close to a hard planar wall, or Ps(r) close to a fixed particle, do display small deviations 
between the different versions of the Rosenfeld functional. Moreover, when calculating the 
depletion potential for size ratios of s = 0.1 or smaller, these deviations are amplified and 
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one observes slightly smaller amplitudes of oscillation for the more sophisticated versions of 
the theory. An example is given in Subsection |V lj| (see, c.f., Fig. O). 



The one-body direct correlation function, defined within density functional theory by 
Eq. (^, can be written as 

4"(r; M) = - E/ - r) (27) 

for the Rosenfeld functional. In the limit where all species have the same radius it is easy to 
check that the weighted densities ria in Eq. and hence $, reduce to the corresponding 
quantities for the pure fluid and, since the weight functions in Eq. ( P7[ ) reduce to the 
weight functions of the pure system, cf^^ reduces to the one-body direct correlation function 
of the pure (s) fluid. The depletion potential is then given by pW{r) = — ln(ps(r)/ps(oo)), 
which is the correct result 0. 
Defining functions as 

the grand potential difference in Eq. (|^) can be written as a sum of convolutions of these 
functions with the weight functions of species b: 

pWtir) = J2[ dV'*"(r') co^ir' - r). (29) 

This expression is valid for arbitrary densities. The dilute limit of species b can now be 
taken, within the Rosenfeld functional, by considering the weighted densities Eq. (pBf ) which 
in this case reduce to 



dilute 



» = E / d'r'p,ir')ujnr-r'), (30) 

where the set of density profiles {pi{r)} that enter Eq. ( pOD is that of the u — 1 component 
fluid, i.e., the one obtained after taking the limit. It follows that the Helmholtz free energy 
in Eq. ( |25| ) and, consequently, the functions \I'" in Eq. ( pSf ) are those of a z/ — 1 component 
mixture. Species b enters into the calculation of the depletion potential, i.e., the dilute limit 
of Eq. (p9D , only through its geometry, i.e., via the weight functions uj^. 
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This feature of the theory becomes especially important if the number of components v is 
small. In the particular case of a binary mixture, = 2, the minimization of the functional 
in the dilute limit reduces to the minimization of the functional of a pure fluid and the 
weighted densities depend only on the density profile Ps(i") of the small spheres: 

nf"*^(r) = I dVp,(r'X(r-rO. (31) 

Although the direct approach of calculating the depletion potential via evaluating grand 
potential differences by brute force requires only a functional describing the pure fluid, the 
above considerations demonstrate that our present approach based on the one-body direct 
correlation function of the big spheres in a sea of small ones requires a functional that 
describes the binary mixture. 



III. RESULTS FROM THE DFT APPROACH AND ASSESSMENT OF THEIR 

ACCURACY 

In this section we examine the accuracy of some of the approximations inherent in the 
present DFT approach by comparing our DFT results for the depletion potential with those 
of simulations and with the predictions of the general asymptotic theory given in Subsec- 
tion 111 C . 



A. Consistency check 

We consider first the results of two separate routes to obtaining the dilute limit for the 
case of a binary hard-sphere mixture. In the first route both components of the mixture are 
treated on equal footing so that one calculates both pb(r) and Ps(i") and obtains Wtij) using 
Eq. (p9D. By requiring the chemical potential of the big spheres p^ to become more and 
more negative, the bulk density pb(oo) of this component approaches zero and the dilute 
limit is taken numerically. For all the mixtures we investigated, a bulk packing fraction of 
the big spheres of r^^ = 10~^ was sufficiently small to ensure that the density profile of the 
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small spheres is indistinguishable from that of a pure fluid at the same rjs. Moreover the 
convergence of Ps(r) and Wt{r) to their limiting values is rather fast; an explicit example 
is given in Fig. 1 of an earlier Letter |T^. Using the second route, employing the weighted 
densities of Eq. (pTj), the dilute limit is taken directly in the functional. In Fig. |I] depletion 
potentials corresponding to both routes are shown for a big hard sphere near a planar wall 
and a size ratio s = Rs/Rb = 0.1. The bulk packing fraction of the small spheres is 77^ = 0.3. 
We find excellent agreement between the two sets of results. The same level of agreement is 
found for a wide range of size ratios s and packing fractions rig. From this we conclude that 
the limit can be taken directly in the functional, which makes the calculations significantly 
easier to perform, and all the results for the depletion potential we present subsequently will 
be based on this route. 



B. Comparison with simulation data 

The results presented in Fig. |^ test the self-consistency of the two routes to the dilute 
limit within the given DFT approach. In order to test the accuracy of approximations 
introduced by employing the Rosenfeld functional the results of the present approach are 
compared with those of simulations. Fortunately some independent sets of simulation results 
for depletion potentials are available for both the sphere-sphere and the wall-sphere case. In 
Fig. 1^ the depletion potentials between two big spheres in a sea of small spheres, at a size 
ratio of s = 0.1 and various packing fractions up to rjs = O.Qtt/G ~ 0.314, obtained from the 



molecular dynamics simulations of Ref. |1J] are compared with results of the present DFT 
approach. The agreement between the latter and the simulations is generally very good. At 
the higher packing fractions small deviations can be seen near contact and near the first 



minimum but the agreement is within the error bars of the simulations ||28| which are not 
indicated here. We note that in Ref. the depletion force was the quantity measured in 
the simulations and the depletion potential was calculated by integrating a smoothed force. 
For a higher packing fraction, rjg = O.Tvr/G ~ 0.367, the agreement between the depletion 
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potential obtained in the simulations of Ref. [|T^] and our present result is poorer (not shown 



in Fig. 0), but for this large value of rjs the error bars of the simulations are probably bigger 
than for small values of rjs ■ 



In Ref. the depletion potential for a single big hard sphere near a planar hard wall 
calculated within DFT was compared with the results of two independent sets of simulations 
for a size ratio s = 0.2 and a packing fraction r]s = 0.3. Very good agreement was found. 
In Fig. 1^ we present a comparison of our results with simulation results from Ref. [jl5l for 



size ratios s = 0.2 (a) and s = 0.1 (b), for various packing fractions of the small spheres 
up to rjs = 0.3. The original simulation results did not oscillate around W = which led 



us to follow the procedure described in Ref. [16| and to shift the data by a small constant 



amount in order to match the contact values with those of our DFT result. We note that 
in the simulations of Ref. [|1^ the depletion force was measured and the depletion potential 



was obtained by integrating the force. Since the data for the force are available only for 
h < hmax, the integral depends on the cut-off hmax-^^ surmise that this cut-off dependence 
is responsible for W{h) not oscillating around zero. The agreement between our DFT results 
and those of the shifted simulation data is very good. The differences probably lie within 
the error bars of the simulations, for all packing fractions when s = 0.2, and for t]s = 0.1 
and Tjs = 0.2 when s = 0.1. However, for rjs = 0.3 and s = 0.1 clear deviations remain 
between our results and those of the simulations. In this case the shifted simulation data 
for the depletion potential are close to the DFT results for h < ag - the height and position 
of the first maximum are the same - but, in contrast to the DFT results, the simulation 
data do not oscillate around zero. Clearly some alternative procedure for interpreting the 
simulation data is required. 

C. Density profiles 

In Fig. ^ the number density profiles of a binary hard-sphere mixture near a planar hard 
wall as obtained from DFT are shown for three size ratios. The packing fraction of the small 
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spheres is r]s = 0.3 and that of the big spheres is r]h = 10~^. The latter is sufficiently low 
that the density proffies of the small spheres, ps{z), shown in Fig. ^(a), are practically equal 
to that corresponding to pure small spheres. Therefore, they are indistinguishable for all 
size ratios. Because of the hard-body interaction between the small spheres and the wall the 
density profile Ps{z) exhibits a discontinuous fall to zero at z = Osjl. The density proffies 
of the big spheres for size ratios s = 0.1 (full line), 0.1333 (dotted line) and 0.2 (dashed 
line) are shown in Fig. §(b). These density proffies do differ significantly for different size 
ratios. The hard-body interaction between the wall and the big spheres does not allow their 
centers to encroach closer than z = o\)j1 and we find that the contact value is very different 
in all three cases (see caption to Fig. ^). We note that the wavelength of the oscillations 
in both Psi^z) and Phi^z) is approximately Og. In order to display the asymptotic behavior 
of these density profiles the logarithm of the difference between each density and its bulk 
value is shown in Fig. |^. For zjos ^ 2 these plots conform very closely to the asymptotic 
form given by Eq. (0). Straight lines joining the maxima have a common slope and the 
distance between adjacent maxima is the same in all cases. Only the amplitudes of the 
oscillations in Ph{z) differ for different values of s. It follows that the decay length, Og ^, and 
the wavelength of the oscillation, 211 /ai, are the same for both density profiles, i.e., for the 
big and the small spheres, and are independent of the size ratio. We have con&med that 
the same values for and ai are obtained from plots of the density profiles of the same 
binary mixture in the presence of a fixed big hard sphere, i.e., our results are consistent 
with Eq. (0). At high packing fractions of the small spheres, e.g., rjg = 0.42, we can easily 
resolve up to 25 damped oscillations. At long range we find the calculated density proffies 
to be in excellent agreement with the predictions of the theory of asymptotic decay. As the 
amplitude of the 24th oscillation is smaller than the amplitude of the ffist one by a factor of 
approximately 5 x 10"^ this attests further to the high numerical accuracy of our results. In 
addition we con&med numerically that the modifications of the Rosenfeld functional which 
we employed lead to the same asymptotic behavior of W{h) as the original functional P? 
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D. Asymptotic behavior 



The asymptotic behavior of the depletion potential calculated within DFT is shown in 
Fig. 1^. For zjos > 2 our results conform very closely to Eq. (^4|): although the amplitude 
of the oscillations depends on s, Wi^z) is characterized by the same, common decay length 
and wavelength l^ajax which describe the density profiles of the mixture. The results 
displayed in Figs. |^ and ^ indicate that the asymptotic behavior of the density profiles 
and of the depletion potential set in at rather small distances from the wall. For wall- 
sphere surface separations of typically z ~ la^^ or even smaller, the asymptotic formulae are 
already remarkably accurate. This is in keeping with the results of earlier studies of the bulk 



pairwise correlation functions of hard- sphere mixtures p2[) where leading-order asymptotics 
were shown to be accurate down to second- nearest neighbor separations. We shall make use 
of this observation in a later section in which we develop an explicit parameterized form for 
the depletion potential. 

As a final examination of the validity of the asymptotic analysis we calculated values of oq 
and a\ from plots of (the logarithm of) the density profiles of the small spheres Ps(-2) near the 
planar hard wall (see Fig. ^) for a range of values for 77^ and various size ratios from s = 0.5 to 
s = 0.1. In accordance with the above statement the results for and ai do not depend on s 
and are shown in Fig. |^ together with the values obtained using the Percus-Yevick result for 
cfj (r) in the pure fluid to solve Eq. (|20| ) for the poles q = ai+iao; we recall that the Rosenfeld 
hard-sphere functional generates the Percus-Yevick two-body direct correlation functions for 
a bulk mixture |]13|. For small packing fractions ris the oscillations are damped very rapidly, 
i.e., the decay length Og ^ is small, so that the numerical determination of the wavelength 
from a density profile is quite difficult. Nevertheless, the level of agreement between the two 
sets of results is very good, for all values of ris that were considered, confirming that the 
DFT results are consistent with the general predictions for the asymptotic behavior. 

Our approach predicts depletion potentials for both the wall-sphere and the sphere- 
sphere case which are in very good agreement with simulations for distances close to contact 
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and which are consistent with predictions of the general theory of the asymptotic decay of 
correlations in hard-sphere mixtures for distances away from contact. From our comparisons 
we conclude that our approach yields accurate results in the whole range of distances, for 
packing fractions up to (at least) rjs = 0.3 and for size ratios down to (at least) s = 0.1. 
We emphasize that the full structure of the depletion potential, which is correctly described 
by the present approach, is not captured by the Asakura-Oosawa approximation or by a 
truncated virial expansion P,^]. 



E. Large asymmetries 

So far it is not apparent how well our present approach will fare for extreme asymmetries, 
i.e., for s ^ 1. The Rosenfeld functional, which is the density functional we apply for all 
of the calculations of the depletion potentials, is designed to treat a mult i- component hard- 
sphere mixture with arbitrary inhomogeneities. Whilst its accuracy in describing the density 



profiles for a pure fluid |T3| and for binary mixtures at moderate packing fractions and 
moderate size ratios has been confirmed by comparison with simulation results, a highly 
asymmetric binary mixture has not yet been studied systematically using this functional. 
Thus it is not known for which size ratios the results calculated with this functional are 
accurate. We recall that the Percus-Yevick approximation becomes increasingly less accurate 
for bulk properties as s — > 0, but here we are interested, in particular, with the reliability 
of our approach for determining depletion potentials. The latter are obtained from density 
profiles, having taken the dilute limit of one of the species [see Eq. (^]. 

In this context it is instructive to consider the depletion potential between a hard sphere 
of radius Ri and one of radius R2 in a sea of small hard spheres of radius Rg at a packing 
fraction rig. This system is formally a mixture of three components in which two are dilute. 
The radius ratio Rg/Ri is chosen such that on the basis of our previous results we know that 
the Rosenfeld functional can treat a mixture of species 1 and s accurately. On the other 
hand, the radius R2 is chosen to be much bigger than Rg and Ri. The depletion potential 
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can be calculated in two different ways. In the first route, sphere 2 (with large radius R2) 
is fixed and enters into the calculation as an external potential for sphere 1 and species s. 
The density profile of the small spheres in the presence of this external potential can be 
calculated, and from it the depletion potential, using the theoretical approach described in 
Sec. |T[ Thus, in this calculation the Rosenfeld functional treats a mixture with a moderate 
size ratio Rg/Ri exposed to an external potential. Therefore we expect these results to be 
very accurate. In the second route, the roles of spheres 1 and 2 are exchanged. The sphere of 
medium radius Ri is fixed and acts as an external potential for the very large sphere 2 and 
the small species s. Now the Rosenfeld functional must treat a very asymmetric mixture. Of 
course, in an exact treatment of this problem it does not matter which sphere is fixed first as 
the depletion potential is simply the difference in the grand potential between a configuration 
in which spheres 1 and 2 are fixed and positioned close to each other and one in which both 
spheres are at infinite separation. The result of an exact treatment cannot depend upon 
which sphere is regarded as an external potential. However, it is not immediately obvious 
that the underlying symmetry is respected in an approximate DFT treatment. At first sight 
one way of calculation might appear to be less demanding on the theory than the other. 

In Fig. 1^ we show the depletion potential between the big spheres when sphere 2 is fixed 
(solid line) and then with sphere 1 fixed (symbols), for a packing fraction of the small spheres 
rjs = 0.3, and Ri = bRg. Two different values of R2 are considered, namely Rs/R2 = 0.02 
(a) and Rs/Ri = 0.01 (b). We find excellent agreement between the results of the two routes 
for both (a) and (b). Only very small differences between the curves can be ascertained and 
these occur for separations close to contact where numerics are most difficult. Note that 
the results in Figs, ^(a) and ^(b) lie very close to each other. This can be understood easily 
when i?2 is the fixed sphere. For cases (a) and (b), R2 ^ Ri = 5Rs and one is effectively 
in the planar- wall limit so that both sets of results lie close to those in Fig. ^(a), with 
rjs = 0.3. From the results shown in Fig. ^ and further comparisons for other values of rjs it 
is evident that the Rosenfeld functional does maintain the required symmetry between 1 and 
2. It is important to understand this. In the low density limit, i.e., if the ternary mixture 
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is considered with the density of the small spheres also approaching zero, the depletion 
potential can be expressed in terms of Mayer / functions and the equivalence of the two 
routes can be verified directly. Starting from the functional given in Eq. (^) we follow the 



for that with sphere 1 fixed. Here fis and /2s are the Mayer / functions between a small 
sphere and sphere 1 and 2, respectively, and it is evident that W2i{r) = Wuir). The 
Rosenfeld functional will reproduce this result for packing fractions 77^ — 0, since it reduces 
to Eq. (^) in this limit. For arbitrary values of 1]^ it is necessary to reconsider the genesis 
of the functional and recognize that although the hard-sphere pairwise potentials ^ij{r) 
between species i and j do not enter explicitly, the functional does respect the equivalence 
of ^ij{r) and $jj(r); the Mayer functions and the weight functions which were used in 
constructing the functional are symmetric w.r.t. i and j. It is straightforward to show 
that the equivalence of W2i{r) and Wi2{r) is guaranteed provided the functional respects 
this symmetry. Thus, the two sets of results shown in Fig. ^a) should agree with each 
other, as should those shown in Fig. H(b). That there are small discrepancies refiects only 
numerical inaccuracies rather than any fundamental shortcoming of the DFT approach. It 
is pleasing that what appear to be two distinct ways of calculating the depletion potential 
yield the same results, even for high degrees of asymmetry. Whether other functionals, not 
based on fundamental measure theory, will respect the symmetry requirements remains to 
be ascertained. Although the present calculations should be regarded as a further test of the 
internal consistency of our approach rather than a formal demonstration that it is accurate 
for extreme asymmetries, the results, when coupled with the excellent agreement between 
theory and simulations for s = 0.1, do suggest that the approach should remain accurate 
for smaller size-ratios. 



derivation of Eq. (|TT]) and obtain 




(32) 



for the depletion potential with sphere 2 fixed and 




(33) 
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IV. DERJAGUIN APPROXIMATION 



In the well-known Derjaguin approximation [3^ the force between two large convex 
bodies is expressed in terms of the interaction energy of two parallel plates. This approximate 
mapping is valid in the limit where the minimal separation of surfaces h is much smaller 
than the radii of curvature and was developed assuming the force between the surfaces can 
be calculated by integration over all interactions between pairs of points of the two bodies. 
Recently the Derjaguin approximation was implemented for the depletion force between two 
big hard spheres in a sea of small hard spheres, employing a truncated virial expansion to 
calculate the excess pressure of the small spheres between planar hard walls, and results were 
compared with simulation data for a size ratio s = 0.1 There is, however, an important 
conceptual difference from earlier applications of the Derjaguin approximation as depletion 
effects are global effects arising from packing of the small spheres and it is not obvious that 
the original derivation remains applicable or what the regime of validity of the approximation 
should be. Some of its limitations were discussed in Ref. where it was argued that the 
Derjaguin approximation should not be reliable for s = 0.1 if the packing fraction rjg > 0.3. 
Here we examine some of the key predictions of the Derjaguin approximation by making 
comparison with results of our DFT approach. From the arguments of Subsec. [111E| it is 
safe to assume that the present DFT approach remains reliable for rather large size ratios 
where the Derjaguin approximation might be expected to be valid. 

There is an elegant scaling relation connecting the depletion force, F{h) = —dW{h)/dh, 
between two big spheres, Fbb{h), in a sea of small spheres with that between a single big 
sphere and a planar hard wall, F^i,{h). In the limit of infinite asymmetry, s —>■ 0, the forces 
are equal except for a factor of 2, i.e., 

2Fbb{h) = F^b{h), (34) 

with h the minimal separation of the surfaces of the two big objects. This scaling relation 
follows directly from the Derjaguin approximation and if it is found to be obeyed it is 
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sometimes inferred |1T0|JTT[] that the Derjaguin approximation itself is vahd. However, it 
was shown that this scahng relation follows from geometrical considerations without 
introducing the explicit Derjaguin approximation. This can be illustrated by comparing the 
explicit Asakura-Oosawa depletion potentials [see Eqs. ( [T^ ) and (|13|)]. In the limit Rb ^ Rs 
both formulae reduce to — |ps(oo)7ri?f,(2i?s — h^, for h < 2Rs, where e = 1 corresponds to 
the sphere-sphere and e = 2 to the wall-sphere case. Thus, achieving the correct scaling 
property in Eq. (|3^) does not prove that the Derjaguin approximation, 

oo 

pDerj^j^^ = -STTiRf, + R,) J dL f,{L) (35) 

h 

is accurate. Here fs{L) is the solvation force, or the excess pressure, for the small-sphere 
fluid confined between two planar parallel hard walls separated by a distance L 

From our DFT calculations we find that the scaling relation is already well-obeyed at 
moderate size ratios. In Fig. | the scaled depletion force = 2i3Fbf,{h)Rl/{Rb + Rs) 

between two big hard spheres (solid line) and that between a single big hard sphere and a 
planar hard wall (□), Pf^^W = [3Fwb{h)R?g/ {R^ + Rg), in a sea of small hard spheres at a 
packing fraction of rjs = 0.3 is shown for size ratios s = 0.1 and 0.02. While small deviations 
from the scaling relation in Eq. (0) are visible close to contact for s = 0.1, these deviations 
have almost disappeared for s = 0.05 (not shown in the figure) and near perfect agreement is 
found for s = 0.02. Note also that the scaled depletion forces corresponding to the different 
values of s lie close to each other. 

An explicit result of the Derjaguin approximation [Eq. (^)] is that the depletion force 
between two big spheres or between a big sphere and a planar wall can be written as p 

F''''^{h) = en{Rb + Rs)ipiVs)ih-2Rs)-^iVs)), h < 2Rs, (36) 

where p{t]s) is the bulk pressure of the small spheres and 7(77^) is twice the surface tension 
of the small-sphere fluid at a planar hard wall. The geometrical factor e is the same as 
in Eq. (|35D. Thus, for a given size ratio the slope of the depletion force predicted by the 
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Derjaguin approximation is constant for h < 2Rs and depends only on the equation of state 
of the small spheres p{r]s)- For the particular case of hard spheres we obtain: 

Oh -<^^ + ^^)4Rl (l-vs)' ' " ^ ^ 

where the quasi-exact Carnahan-Starling equation of state |31|] was used. 

However, the depletion forces calculated within the present approach show a qualitatively 
different behavior from that predicted by Eq. (|37|) . It was found that even for small size 
ratios (s < 0.05), only in the limit r]s 0, in which the Asakura-Oosawa approximation 
becomes exact, there is agreement between the Derjaguin approximation and the results of 
our approach. The depletion force calculated at a packing fraction of rjs = 0.3 does not 
have constant slope for h < 2Rs (see Fig. This is in clear contradiction to Eq. (p7|). 
Simulation results for the depletion force also exhibit non-constant slopes for h < 2Rs. 



Another prediction of the Derjaguin approximation in Eq. (|36D is that the contact value 
W{0) of the depletion potential can be expressed simply in terms of the equation of state 
p{i]s) and the surface tension 7(^7^)- Using the Carnahan-Starling result for p{ris) |^T| and 
the scaled particle result for 7(^75) [0 we obtain B: 



= ms ^''^ 

which becomes positive at high packing fractions of the small spheres 0. Provided rjs < 0.2 
the contact values from Eq. (|38D are in reasonable agreement with the results of our DFT 
approach for small values of s. However, the contact values obtained from Eq. (|38[) change 
sign at rjs ~ 0.3532, which is in complete contradiction to the results of the present approach 
where we find negative contact values for all packing fractions and all size ratios s under 
consideration p3| . 

In Ref. [0 it was shown that a third-order virial expansion (in powers of rjs) for the 
depletion potential calculated within the Derjaguin approximation does not yield positive 
contact values. However, expansion to fourth or fifth-order shows a qualitatively different 
behavior from third order and already indicates the onset of positive Vr(0). Thus, in keep- 
ing with Ref. we conclude that the Derjaguin approximation is not very useful for the 
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calculation of depletion forces. The good level of agreement, observed for h < 2Rs, between 
the results of the third-order virial expansion ||^ and those of simulation [|T^ for s = 0.1 
should be regarded fortuitous. 



V. APPLICATIONS 



A. A parametrized form for the depletion potential of hard spheres 



As mentioned in Subsec. |1I A| , recent studies of correlation functions and phase equilibria 



of highly asymmetric binary mixtures have shown that it is very advantageous to map 
such mixtures onto effective one-component fluids [^0. The effective pairwise potential 
between the big particles is then the bare pair potential between two big particles plus the 



depletion potential [see Eq. (14)]. Thus, in calculating the phase behavior of binary hard- 
sphere mixtures it is necessary to adopt a specific form for the depletion potential between 
two big hard spheres. Previous simulation studies of binary mixtures have employed 
the simplified third-order virial expansion formula given by Gotzelmann et al. and the 
same potential has been used in a perturbation theory treatment of the phase behavior 
0]. Although this formula is convenient for global investigations of phase behavior, as the 
depletion potential is given explicitly as a function of rjs, clearly it would be valuable to 
have a simple, parametrized form for the depletion potential that (i) is better founded than 
the formula provided by Gotzelmann et al. and (ii) captures the correct intermediate and 
long-range oscillatory structure as well as the important short-range features. Note that in 
Refs. 1^ and |^ the effective pair potential was set equal to zero for separations h > 2Rs. 

We have used depletion potentials calculated within the present DFT approach, for a 
single big sphere near a planar hard wall and for two big spheres, to develop a suitable 
parameterization scheme. Although this parameterization is fairly simple it yields rather 
accurate fits. The depletion potential close to contact is fitted by a polynomial and is 
continued by the known asymptotic behavior. In the following the variable x measures the 
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minimal distance from contact in units of the small sphere diameter as-, i.e., x = h/ag. These 
parametrized depletion potentials W are also scaled: the actual potentials W are recovered 
by multiplying by a factor of e{Rb + Rs) / (2i?s) with e = 2 for the wall-sphere and e = 1 for 
the sphere-sphere potential: 

W = <^^±Mw. (39) 

Between contact at a; = and the location xq of the first maximum the scaled depletion 
potential is fitted by a cubic polynomial: 

PW{x,r]s) = a{r]s) +b{r]s) X + c{r]s) x^ + d{r]s) x^, x < Xq, (40) 

where the coefficients a, b, c and d are functions of the packing fraction of the small spheres 
?7s. More details of this polynomial and the determination of the coefficients are presented 
in the appendix. 

In order to obtain the depletion potential for a; > Xq we assume that the asymptotic 
decay already sets in at the point xq. This assumption is supported by the results presented 
in Fig. ^. Thus, for x > Xq we adopt the form [c.f. Eq. 



P^asympti.^^ Vs) = A^iVs) exp(-ao(?7s)cT,x) COs(ai (r7,)cr,X - QujiVs)), X > Xq, (41) 

for the scaled depletion potential between a wall and a sphere and [c.f. Eq. (l2l|) ] 

A (rjs) 

PW^sympt{x,Vs) = _^ ' exp{-ao{r]s)cTsx) cos{ai{r]s)cTsX - Qp{r]s)), x > Xq, (42) 
s + X 

for the potential between two spheres. The denominator in Eq. (|42|) measures the separation 
af, + h between the centers of the spheres in units of ag- Both forms contain the functions 
aQ{ris) and ai{r]s), which can be calculated from the Percus-Yevick bulk pair direct correlation 
function cis^(r) (see Subsec. |111 D| and Fig. ^. The amplitudes Aj{ris) and phases 9j(?7s), 



j = p,w, are chosen so that the depletion potential and its first derivative are continuous at 
Xq. Ap{r]s) and Qp{ris) are weakly dependent on the size ratio s. 

With this prescription the scaled depletion potential is completely determined. For 



a given packing fraction rjs the coefficients a, 6, c and d are given by Eq. ( |A7| ) and the 
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position of the first maximum can be calculated from Eq. Using those values as input, 

the amplitude A and the phase G of the asymptotic decay are readily obtained from either 
Eqs. and ( [A3| ) or Eqs. ( [A6| ) and Thus in this parametrization the scaled depletion 

potential has the form 

{a + b X + c x'^ + d x^, x < Xn 
(43) 
P^asympti^yVs): X > Xq. 

In Fig. 0(a) fits (lines) of the form given by Eq. (^) are compared with the scaled 
depletion potentials between a big hard sphere and a hard wall calculated within DFT 
(symbols) for a size ratio s = 0.1. Although the fit is relatively simple, its accuracy is high. 
The position of the first maximum, which depends sensitively on the packing fraction rjs, is 
reproduced very accurately. The value PWq = pW{xo) of the potential at the first maximum 
is also given quite accurately, and only for r]s = 0.3 are small deviations of the fit from the 
full DFT results visible. Clearly the full structure of the depletion potential is reproduced 
well by this parametrization. In order to demonstrate the wide range of applicability of this 
parametrization in Fig. |lO|(b) we show a comparison of the parametrized scaled depletion 
potential (□) for a packing fraction ris = 0.3 with scaled DFT results (lines) for the depletion 
potential between two spheres and size ratios s = 0.2 and s = 0.05. Although for these size 
ratios the scaling relation Eq. (Q) is not satisfied particularly accurately, the agreement 
between our parametrization and the DFT results is rather good. This gives us confidence 
that we have developed a satisfactory parametrized form for the depletion potential which 
properly incorporates all essential features. 



B. Oscillatory depletion potential at high packing fractions of the small spheres 

In a recent experiment by Crocker et al. |^ the equilibrium probability distribution p{r) 
for two (big) PMMA (polymethylmethacrylate) spheres of diameter cxb = l.lfim immersed 
in a sea of (small) polystyrene spheres of diameter as = 83 nm was measured using line- 
scanned optical tweezers and digital videomicroscopy at various packing fractions in the 
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range between rjs = 0.04 and ris = 0.42. The solvent contains added salt and surfactant 
to prevent colloidal aggregation and the 'bare' interactions between the colloidal particles 
are expected to be screened Coulombic repulsion with a screening length of about 3 nm 
0. Since the latter is small compared with the colloid diameters the bare interactions 
can be regarded, to good approximation, as hard-sphere-like. The depletion potentials 
f3W{r) = —ln{p{r)/p{oo)) (see Subsec. J| A) obtained from these experiments are shown in 
Fig. 1 of Ref. 0. At low packing fractions, rjs = 0.04 and 0.07, rather good agreement with 
the results of the Asakura-Oosawa approximation was found, after taking into account the 
effects of limited spatial resolution of the optical instruments. For ris = 0.15 and 0.21 the 
measured depletion potential displayed a pronounced repulsive barrier. For higher packing 
fractions, i.e., r]s = 0.26, 0.34 and 0.42 damped oscillations were observed, these being 
particularly pronounced for the two highest packing fractions for which three maxima are 
clearly visible. Reference appears to be the first report of an experimental observation 
of an oscillatory depletion potential and, indeed, of a repulsive contribution arising from 
purely entropic or packing effects 



Motivated by these experiments we consider a binary hard-sphere mixture in the dilute 
limit with a size ratio s = 0.0755, as in the experiment. (We do not attempt to include the 
increase of the effective radius of the spheres arising from screened Coulomb repulsion and, in 
keeping with the authors of Ref. 0, we do not include any dispersion forces.) As previously, 
the depletion potential between two big spheres is calculated using Eq. (^) in the dilute 
limit. The functions are functionals of Ps(r), the density profiles of the small spheres close 
to a big sphere fixed at the origin, which depend only on the radial distance r. The results 
are shown in Fig. ^ for the same values of ris as in the experiments. It is encouraging to 
find that the theoretical and experimental results have many common features. As expected, 
the calculated oscillations become much more pronounced as ris increases. The wavelength 
decreases slowly and the decay length of the envelope increases rapidly with ris - as predicted 
by the theory of asymptotic decay (see Fig. |^). The experimental data are consistent with 
both observations. Moreover the wavelength of the oscillations for ris = 0.34 is close to 
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(Ts = 83 nm in theory and experiment. For r]s = 0.42 both theory and experiment yield a 
shghtly smaller wavelength. The amplitude of the calculated oscillations is larger than in the 
experiment. However, we emphasize that we made no attempt to take into account effects of 
instrumental resolution or the polydispersity of the small polystyrene particles. Nor have we 
attempted to include the effects of the softness of the inter-particle potentials and any non- 
additivity of the effective diameters; both are likely to lead to a reduction in the amplitude 
of the oscillations. The qualitative agreement between the experimental results and those 
of our calculations persuades us that the hard-sphere model is an appropriate starting point 
for describing the colloidal system and that the observed oscillations do reflect the packing 
of the small spheres - as inferred in Ref. 0. 

Signiflcant deviations between our results and the experimental ones do occur, at large 
rjs, for separations near contact or near the flrst maximum in the depletion potential. Our 
results imply that the height of the flrst maximum and the magnitude of the contact value 
I W{0) I are larger than the experimental ones by about a factor of two for rjs = 0.34. Although 
the source of these differences may well reside in the experimental situation it is important 
to check that the particular DFT which we employ is performing reliably at these high values 
of ris. It is precisely this regime of high density and very strong conflnement of the small 
spheres where differences between the various DFT theories, i.e., the improvements on the 
original Rosenfeld version, might reveal themselves. These circumstances are reminiscent 



of those investigated by Gonzalez et al. |3^ in their DFT studies of hard spheres in small 
spherical cavities. Those authors were able to ascertain that the improved theories fared 
better than the original version under conditions of extreme conflnement. 

To this end we repeated our calculations of the depletion potential with the improved 
versions of the Rosenfeld functional that can account for the freezing transition [^]. At 
packing fractions 7]^ < 0.3 we obtained, as stated earlier, results almost identical to those 
of the original functional. At higher packing fractions, however, we flnd that the amplitude 
of the oscillations is slightly smaller than those obtained from the original functional. This 
is illustrated in Fig. for rjs = 0.42, using the interpolation form of the functional ||27 
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(dotted line). The antisymmetrized version of the functional, with g = 3 yields a 
depletion potential very close to that of the interpolation form. In view of the smallness 
of these deviations the discrepancies between the experimental findings and the theoretical 
results at high rjs cannot be blamed on the performance of the DFT but most probably 
reside in differences between the actual experimental sample and the model of hard spheres. 

VI. SUMMARY AND DISCUSSION 

In this paper we have developed a versatile theory for determining the depletion potential 
in general fluid mixtures. Our approach requires only the knowledge of the equilibrium 
density profile Ps(r) of the small particles before the big (test) particle is inserted, i.e., Ps(r) 
has the symmetry of the external potential. If the latter is exerted by a fixed particle or by 
a planar wall then in these cases Ps(r) simplifies to functions Ps(r) or Ps{z) of one variable. 
Since a one-dimensional profile can be calculated very accurately, the resulting depletion 
potentials can be obtained without the numerical complications and limitations that are 



inherent in brute-force DFT [^. The latter requires the calculation of the local density of 
the small particles around the big particles in the presence of the external potential ^Tf or 
the calculation of the total free energy as a function of the separation of the big particles 



37| ; both calculations require considerable numerical effort due to the reduced symmetry 
of the density distributions. We have employed our approach in a comprehensive study 
of the depletion potential for hard-sphere systems, using Rosenfeld's fundamental measure 
functional. The main conclusions which emerge from our study are as follows: 

1. The depletion potential can be obtained by considering a liquid mixture in the limit 
of vanishing concentration of one of the species. Two different ways to implement this 
limit lead to the same result (Fig. 

2. Detailed comparison of our results with those of simulations, for both sphere-sphere 
and (planar) wall-sphere depletion potentials (see Figs.^ and ^, demonstrate that 
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the theory is very accurate for size ratios s = Rs/Rb as small as 0.1 and for packing 
fractions rjs as large as 0.3. These are the most extreme cases for which reliable 
simulation data are presently available. The theory describes accurately the short- 
ranged depletion attraction, the first repulsive barrier and the subsequent oscillations 
in the depletion potential. 

3. By performing consistency checks we argue that at least up to moderate packing 
fractions the predictions of the Rosenfeld DFT for depletion should be quantitatively 
reliable even for large asymmetries between the sizes of the solvent and the solute 
particles (Fig. |). Subsection |III E| provides a theoretical understanding of this feature 



of our DFT approach. 

4. Extensions of the Rosenfeld functional ||2^ yield very similar results (see Fig. |lT]) for 
the cases we have studied. It would be of considerable interest to test the performance 
of the proposed functionals against simulation data for smaller size ratios and for 
higher values of ris, for which more extreme packing constraints might discriminate 
between the various functionals. 

5. Our DFT approach incorporates the correct, exponentially damped, oscillatory asymp- 
totic {h oo) decay of the depletion potential W{h). This is inherent in the con- 
struction of the theory, is preserved by the approximate Rosenfeld functional and is 
exhibited explicitly by the numerical results (Fig. The decay length a^^ of the 
oscillations increases and the wavelength 27r/ai decreases with increasing ris (Fig. ^ 
but these quantities are independent of the size ratio s. The same values for ao and 
ai characterize the oscillatory decay towards the bulk values of the number density 
profiles of hard-sphere mixtures near a hard wall when the packing fraction of the big 
spheres is very small (Figs. ^ and 

6. We have developed simple parametrization schemes for the depletion potential between 
big hard sphere and a planar wall and that between two big hard spheres which provide 
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accurate fits to our DFT results (see Fig. |^). The fitting procedure makes use of the 
fact that leading asymptotic behavior of W{h oo) provides an accurate account of 
the oscillatory structure of the depletion potential at intermediate separations as well 
as at longest range. Such parametrizations are designed to provide a more accurate 
alternative to the third-order virial expansion formula given by Gotzelmann et al. 

Since these parametrization can be easily implemented we recommend that they 
should be employed in subsequent studies of the phase behavior of highly- asymmetric 
binary hard-sphere mixtures of the type reported in Refs. and 



7. In Sec. |^ we investigated the regime of validity of the Derjaguin approximation 
[Eq. (^)] for the depletion potential and showed that this fails, for all but the smallest 
packing fractions rjs, for which the depletion potential reduces to the Asakura-Oosawa 
result. However, the scaling relation Eq. (|34D connecting the depletion force between 
two big spheres to that between a big sphere and a planar wall - which is predicted by 
the Derjaguin approximation but which also follows from geometrical considerations 
- does remain valid even at moderate size ratios (Fig. 

We conclude with several remarks concerning the accuracy and usefulness of our approach. 
One might be surprised that a DFT which corresponds to the Percus-Yevick theory for the 
bulk mixture (the Rosenfeld functional yields the same bulk free-energy density and bulk 
pair direct correlation function) performs so well for small size ratios, for which it is known 
that Percus-Yevick theory becomes inaccurate. For example, it fails to predict the fiuid-fiuid 
spinodals for hard-sphere mixtures. However, our present approach involves only the calcu- 
lation of a one-body direct correlation function c[^^(r; {/ij}) and, therefore, the determination 
of one-body density profiles. The minimization of approximate functionals can yield rather 
accurate one-body profiles in spite of the limitations of the underlying approximations; e.g., 
this is the reason why the test particle route to the bulk radial distribution function g{r) 
is very successful within DFT [p!3| , |38[] . Furthermore, in determining the depletion potential 
W{r) we require only solutions of the Euler-Lagrange equation for Ps(r) in the limit where 
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Ph — > 0, i.e., in the absence of the big particles. The DFT is hkely to be more accurate in this 
hmiting regime than for a mixture concentrated in all species. We emphasize that taking the 
dilute limit of the big particles numerically, i.e., working at non-zero but very small values 
of r/b, involves more computation than taking the limit directly in the functional. Moreover, 
caution should be exercised in hard-sphere mixtures with extreme size ratios, s < 0.1, at 
high packing fractions ris of the small spheres since the fluid-solid phase boundary already 
occurs at very low packing fractions T]b of the big spheres [@|. The fluid-solid coexistence 
region is avoided if the dilute limit is taken directly []39| . 

Our procedure for calculating Wt{r) at arbitrary concentrations of the big particles might 
prove useful for interpreting (future) measurements of the effective interaction potential when 



the mixture is not in the dilute limit. Figure 1 of Ref. [|T6[ illustrates how the wall- sphere 
potential Wt{z) varies with the big sphere packing fraction rjf, for a mixture with size ratio 
s = 0.2 and r]s = 0.2. For rjb = 0.025, Wt{z) already differs by a few percent from its dilute 
limit W{z). 

It is possible to calculate the depletion potential by using as input density profiles ob- 
tained by other means. In particular one might take simulation data for Ps(r), computed 
in the absence of the big test particle, and insert these into Eq. (^Tj) to determine the 
weighted densities. Although such a procedure does not offer the appeal of a self-consistent 
approach in which both the equilibrium density profiles and the depletion potential are cal- 
culated within the same framework, in practice this could be a profitable route for complex 
geometries where a direct simulation of the depletion potential or force is very difficult. 

Finally we mention that the techniques we have developed here are not restricted to 
additive, binary hard-sphere mixtures. Our general approach to the calculation of depletion 
potentials can be applied to hard-sphere mixtures with non-additive diameters, to ternary 
mixtures and to systems where the interparticle potentials are soft. 
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APPENDIX A: PARAMETERIZING THE DEPLETION POTENTIAL 

In the range between contact, x = h/as = 0, and the position xq of the first maximum the 
scaled depletion potential W is parametrized by a cubic polynomial [Eq. (|40|)], which is the 
simplest polynomial fit that remains accurate close to contact. Since pW{x = 0, rjs) = a{ris), 
the first coefficient is the contact value of the depletion potential. The position xo and the 
height Wo of the first maximum can be obtained easily by differentiating Eq. (^Ol). The 
cubic polynomial has two extrema, with the maximum located at 



, , c + \/c^ - 3bd , . ^ , 

MVs) = ^ (Al) 



and a maximal value of 



mMs) = PWix = xoiVs),Vs) = —■ (A2) 

Beyond the position of the first maximum of the depletion potential the parametrized 
form is continued by imposing the known asymptotic behavior for large h. The asymptotic 
behaviors of the wall-sphere and sphere-sphere depletion potential are slightly different and 
must be considered separately. For the wall-sphere depletion potential the asymptotic be- 
havior is given by Eq. (^TJ) and the amplitude and the phase 9^ are chosen such that the 
function and its first derivative are continuous at Xq. From the requirement of a continuous 



XUIVC CLlj IjllC lllOL llid-A-llllUlli, 1.^., 

to be 



derivative at the first maximum, i.e., 



= 0, the phase can be determined 

x=xo 
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Qwiris) = aiasXo + arccos ( ) . (A3) 

Vag + aij 

From the requirement that the function is continuous at Xq, i.e., W^^y^p^{xo,ris) = Wo{ris), 
together with the phase from Eq. (|A3|) , the amphtude of the asymptotic decay follows as 



AUVs) = pWoexp{aoasXo)\l ^^\"'\ (A4) 



A similar calculation for the sphere-sphere case using Eq. (g^) leads to slightly different 
expressions for the phase, 

QM = a.a^xo + arccos ( aA^.S + a,) \ ^^^^ 

+ 2ao(xoa, + at) + (a^ + af)(xoo-s + ) 

and the amplitude 



A,{r^s) = ^^^ ^^P (^o^-^o) L ^ 2ao{xoa, + a,) + (ag + af)(xoa, + a,)^. (A6) 
aiCTs V 

Unlike 9^^, and A^, Qp and Ap depend (weakly) on the size ratio s = ag/cTb- 

The coefficients a, b, c and d are fitted to depletion potentials calculated within DPT. 
Scaled depletion potentials W obtained for a big hard sphere near a planar hard wall, for a 
size ratio s = 0.1, are used in the range < r^s < 0.3. The resulting coefficients are given by 

a{7]s) = -2.909 r]s, 

biris) = 6.916 77, - 4.616 77^ + 78.856 rit 

(A7) 

c{r],) = -4.512 r]s + 15.860 r^f _ 93.224 r]^ 

d{ris) = -ris exp(-1.734 + 8.957 77, + 1.595 t/^). 
There is no particular significance in the chosen form of parametrization but we note that 
the contact values of the scaled depletion potential l3W{0,ris) = a{r]s) are linear in 77, for 
this choice of parametrization. It is interesting that the coefficient —2.909 is rather close 
to the value —3 obtained from the Asakura-Oosawa result (valid as rjg —>■ 0) in the limit of 
small size-ratios, see Eqs. (|T2D and (O) and also Eq. (pSj). 



The quantities and ai are obtained by solving Eq. (pOD, using the Percus-Yevick 
pair direct correlation function c^^J(r). The results are shown in Pig. |^ In the range 
0.05 < ris < 0.4 they can be fitted accurately by 
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ao{r]s) (Ts = 4.674 exp(-3.935 77^) + 3.536 exp(-56.270 77^) (A8) 

and 

ai{rjs) as = -0.682 exp(-24.697 77^) + 4.720 + 4.450 77^. (A9) 

These formulae specify all the ingredients for determining the parametrized form of the 
depletion potential. 
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FIGURES 




h/03 

FIG. 1. Comparison of the depletion potential for a big hard sphere near a planar hard wall, 
calculated via two different routes. These correspond to taking the dilute limit directly in the 
functional (solid line), and to taking the dilute limit numerically with a packing fraction of the big 
spheres of rjb = (□). In both cases the size ratio is s = Rg/Rb = 0.1 and the packing fraction 
of the small hard spheres is rjs = 0.3. h is the separation between the wall and the surface of the 
big sphere; ag = 2Rs is the diameter of the small spheres. 



43 



CO. 



2 ■- 



■■ 



-2 ■■ 








= 0.2 71/6 
= 0.4 71/6 
= 0.5 71/6 
r|s = 0.6 71/6 



H — I — I — I — I — I — I — I — I — I — I — I — I — I — I — |- 



0.5 



1 1.5 

h/o. 



2.5 



FIG. 2. The depletion potential between two big hard spheres in a sea of small hard spheres 
calculated for various packing fractions ijg of the small spheres and a fixed size ratio s = 0.1. We 



compare simulation results (symbols) from Ref. \14\ with results from our DFT approach (solid 
lines), h is the separation between the surfaces of the two big spheres and as is the diameter of 
the small spheres. 
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FIG. 3. The depletion potential between a big hard sphere and a planar hard wall calculated 
for various packing fractions ijg of the small hard spheres and size ratios s = 0.2 (a) and s = 0.1 
(b). We compare processed simulation data from Ref. ||l^ (symbols) with those of our calculations 
(solid lines). The only significant deviations occur in (b) for rjs = 0.3 (see text), h is the separation 
between the wall and the surface of the big sphere; as is the diameter of the small spheres. 
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FIG. 4. The number density profiles of the two components of a binary hard-sphere mixture 

near a planar hard wall as obtained from DFT. The reservoir packing fraction of the small spheres 
is rjs = 0.3 and that of the big spheres rjb = 10~'^. For the three different size ratios s = 0.1, 
0.1333, and 0.2 the profiles of the small spheres (a) are indistinguishable while the profiles of the 

big spheres (b) differ considerably. The contact values of the density profiles of the big spheres are 



Pb{z 



0.0279, 0.2518, and 2.1814 for s=0.2, 0.1333 and 0.1, respectively. z = denotes 



the position of the wall and ag is the diameter of the small spheres. 
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FIG. 5. The asymptotic decay of the number density profiles of the two components of a 
binary hard-sphere mixture near a planar hard wall for three different size ratios s. The packing 
fraction of the small spheres is rjg = 0.3 and that of the big spheres is % = 10~^. In each case the 
natural logarithm of the modulus of the density profile minus the corresponding bulk density is 
plotted. The upper lines denote the density profiles of the small spheres; these are practically the 
same for all s. The lower lines denote the density profiles of the big spheres and it can be seen that 
the amplitude of the oscillations does depend on the size ratio s. However, for z/cg > 2 the same 
characteristic decay length Oq ^ and wavelength of oscillation lir/ai characterize the decay of both 
density profiles; i.e., for the big and the small spheres (see text). Note that the density profiles in 
Fig. 1^ have been shifted so that here z measures the distance from contact, as is the diameter of 
the small spheres. 



47 




— I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 

1 2 3 4 5 

Z/O3 

FIG. 6. The asymptotic decay of the depletion potentials W{z) for a big hard sphere near a 
planar hard wall for the same parameters as shown in Fig. |5|, i.e., r]s = 0.3 and three size ratios 
s. The same decay length ^ and wavelength 27r/ai that determine the asymptotic decay of the 
density profiles determine the decay of W{z) (see text). Only the amplitude of the oscillations 
depend on s. In each case z measures the distance from contact, fx^ is the diameter of the small 
spheres. 
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FIG. 7. Comparison of (a) the inverse decay length gq and (b) the inverse wavelength 
(27r/ai)^^ as determined from the theory of asymptotic decay of the bulk pairwise correlation 
[see Eq. ([20|)] using the Percus-Yevick two body direct correlation function c'^J (r) (solid line) with 
the corresponding results obtained from density profiles calculated using the Rosenfeld functional 
(□) for hard-sphere fluids near a planar hard wall. r]s is the packing fraction and ag is the diameter 
of the small spheres. 
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FIG. 8. The depletion potential between a hard sphere 1 with radius Ri = 5Rs and a hard 
sphere 2 with radius R2 = 50Rs (a) and R2 = lOORg (b), in a sea of small hard spheres with 
radius Rg and packing fraction r)s = 0.3. The solid line denotes the depletion potential calculated 
by fixing sphere 2 first, so that 2 enters into the calculation as an external potential, whereas the 
symbols denote the depletion potential obtained with sphere 1 acting as the external potential. 
The two routes should lead to the same results (see text) . h is the separation between the surfaces 
of the spheres 1 and 2 and ag = 2Rs. 
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FIG. 9. The scaled depletion force between two big hard spheres Pf*{h) = 2PF{h)R'^/{Rb+Rs) 
(solid line) and between a single big hard sphere and a planar hard wall (3f*{h) = 
(3F{h)Rl/{Ri) + Rs) (□) in a sea of small hard spheres at a packing fraction rjs = 0.3 for size 
ratios s = 0.1 (a) and s = 0.02 (b). For s = 0.02, the scaling relation in Eq. ( p^ ) is obeyed almost 
perfectly, h is the separation between the surfaces of the big spheres or between the wall and the 
surface of the big sphere. CTg = 2Rs is the diameter of the small spheres. 
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FIG. 10. (a) Comparison of the scaled wall-sphere depletion potential [Eq. (pOl)] for various 
packing fractions rjs of the small hard spheres and size ratio s = 0.1 as calculated fully within 
DFT (symbols) and as given by the parametrization of Eq. ( ^3| ) (lines), (b) Comparison of the 
scaled sphere-sphere depletion potentials for a packing fraction rjs = 0.3 and size ratios s = 0.2 



and s = 0.05 as calculated within DFT (lines) and as given by the parametrization of Eq. (43) 
(□). Differences between the parametrized results for s = 0.2 and s = 0.05 are not visible, h is the 
separation between the wall and the surface of the big sphere or between the surfaces of the two 
big spheres, as is the diameter of the small spheres. 
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FIG. 11. The depletion potential between two big hard spheres in a sea of small hard spheres 
at various values of the small sphere packing fraction rjs as obtained from the original Rosenfeld 
functional. In order to mimic the experiment of Ref. Q the diameters were chosen to be cjfe = l.l//m 
and as = 0.083;um so that the size ratio is s = 0.0755. h measures the separation between the 
surfaces of the big spheres. Note that for display purposes each curve has been shifted downward 
by a constant amount; W{h) oscillates around zero as /i — > oo. The dotted line for rjg = 0.42 corre- 
sponds to the depletion potential calculated with the modified interpolation form of the Rosenfeld 



functional |27|. 
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